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A quasi-three-dimensional analysis has been developed for unsteady rotor- 
stator interaction in turbomachinery. The analysis solves the unsteady Euler 
or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It 
accounts for the effects of rotation, radius change, and stream-surface thick- 
ness. The Baldwin -Lomax eddy-viscosity model is used for turbulent flows. 

The equations are integrated in time using an explicit four-stage Runge-Kutta 
scheme with a constant time step. Implicit residual smoothing is used to in- 
crease the stability limit of the time-accurate computations. The scheme is 
described, and stability and accuracy analyses are given. 

Results are shown for the stage of the space shuttle main engine high 
pressure fuel turbopump. Two stators and three rotors were used to model the 
41:63 blade ratio of the actual machine. Implicit residual smoothing was used 
to increase the time step limit of the explicit scheme by a factor of six with 
negligible differences in the unsteady results. We feel that the implicitly 
smoothed Runge-Kutta scheme is easily competitive with implicit schemes for 
unsteady flows whi'e retaining the simplicity of an explicit scheme. 
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UNSTEADY ROTOR-STATOR INTERACTION CODE 

DEVELOPED DY P. C. E. JORGENSON k R. V. CIIIMA 


DESCRIPTION 

• UNSTEADY THIN-LAYER NAVIER-STOKES SOLVER 
FOR ROTOR-STATOR INTERACTION 

FEATURES 


• QUASI-3-D FORMULATION 

INCLUDES ROTATION, RADIUS CHANGE, k STREAM SURFACE THICKNESS 

• C-TYPE GRIDS OVERLAP AT INTERFACE 

• MULTI PASSAGE CAPABILITY 

• EXPLICIT RUNGE-KUTTA SCHEME + IMPLICIT RESIDUAL SMOOTHING 
UNCONDITIONALLY STABLE (PRACTICAL LIMIT CFl = 0(20)) 

2nd ORDER TIME ACCURATE 

FIRST USE OF IMPLICIT SMOOTHING FOR UNSTEADY PROBLEMS 
RESULTS 

• SPACE SHUTTLE MAIN ENGINE TURBINE 


CD-89-39438 


QUASI-THREE-DIMENSIONAL STREAM SURFACE 

ROTATION 
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QUASI-3-D THIN LAYER NAVIER-STOKES EQUATIONS 
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WHERE 


wE “ e m v m + f e w 0, 
wn - n m v m + n 0 w 9 


CONTRAVARIANT VELOCITY COMPONENTS IN RELATIVE SYSTEM 


w 0 - v 0 - rfi - RELATIVE TANGENTIAL VELOCITY 

Q - BLADE ROTATION SPEED 

2 -1 -1 
K 2 - (pv 0 + p-Re o 2 2) + (p ■* Re O33) 6 m h(h 

- CENTRIFUGAL FORCE TERM + pd m (AREA) 


- E^/r 

: 0 -METRICS SCALED BY 1 /r 

He - ne /r 

Uj-rh/J : JACOBIAN SCALED BY AREA 
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ROTOR-STATOR INTERACTION CODE DETAILS 

INITIAL CONDITIONS 

ANALYTIC 1-D SOLUTION OF FLOW EQUATIONS WITH AREA CHANGE 

BOUNDARY CONDITIONS 

INLET 

P 0 , Tq, AND v 0 SPECIFIED 
R' = v m - ^ EXTRAPOLATED 

III y-l 

EXIT 

P SPECIFIED 

3 CONSERVATION VARIABLES EXTRAPOLATED 
WALLS 

INVISCID FLOW - TANGENCY, p EXTRAPOLATED 
VISCOUS FLOW ■ NO-SUP, T w SPECIFIED 
NORMAL MOMENTUM EQUATION FOR PRESSURE 

PERIODIC BOUNDARIES 

SOLVED LIKE INTERIOR POINTS 

CD-89-3944 1 
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MULTISTAGE RUNGE-KUTTA ALGORITHM 

GOVERNING EQ UATIONS 
dtq= -J[Rj ~ [Ry + D)} 

/?/ = INVISCID RESIDUAL 
Rv = VISCOUS RESIDUAL 
D = ARTIFICIAL DISSIPATION TERM 
MULTISTAGE SCHEME 


<70 — Qn 


<7i = 9o “ J At [Rj g 0 - [Rv + D) go] 


<lk - q o - <*kJ At \Ri gjt„i - [Ry + 1?) g 0 ] 

<7n+l ~ flJb 

7?v D EVALUATED AT FIRST STAGE ONLY 
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ARTIFICIAL DISSIPATION 

NONCONSERVATIVE VERSION OF JAMESON FORMULATION 
Dq = (D( + D n )q 

f -DIRECTION OPERATOR 

D V 7 = C ~ v *1 ( {{«) 

WHERE: 

C- — - — 

J&tij 

(LOCAL At,,/ MINIMIZES DISSIPATION) 

Vj — /i 2 max (('i-f i , , «'i_ i ) 

IT = max (0, pi — V 2 ) 

AND 

v - ElLIr ~ 2P '-i + p -wl 

|P i+ i, / + 2P f , J +F i _ Ii; | 

/< 2 = 0(1) 

/‘4=0(i) 
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IMPLICIT RESIDUAL SMOOTHING 


USE A TIME STEP GREATER THAN THE STABILITY LIMIT 
MAINTAIN STABILITY BY SMOOTHING THE RESIDUAL IMPLICITLY 

REWRITE STAGE (K) OF MULTISTAGE SCHEME AS 
A q W = q W - q W = -aj&t - D(J W ] 

IMPLICIT SMOOTHING STEP 

(1 - tf«)(l - £c5„„)a7u = A q w 

q( k l = ql°) + AyW 
UNCONDITIONALLY STABLE IF 


WHERE 

CFL‘ IS COURANT LIMIT OF THE UNSMOOTIIED SCHEME 
CFL IS THE LARGER OPERATING COURANT NUMBER 
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IMPLICIT RESIDUAL SMOOTHING 


TWO-DIMENSIONAL LOCAL SMOOTHING PARAMETER 



WHERE 

CFL * IS COURANT LIMIT OF THE UNSMOOTIIED SCHEME 
CFLij IS THE LOCAL COURANT NUMBER 


CD-89-39445 
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IMPLICIT RESIDUAL SMOOTHING 


FEATURES 

• IRS DOES NOT EFFECT FORMAL TIME ACCURACY 
OF R-K SCHEME 

• SMOOTHING ONLY OCCURS WHERE LOCAL COURANT NUMBER CFL 
IS GREATER THAN THE COURANT LIMIT OF THE R-K SCHEME CFL* 
(EG - VISCOUS REGION) 

• CAN USE CFL = 0(10 - 20) IN VISCOUS REGIONS WHILE CFL = 0(A) 
IN INVISCID CORE 

• eh EASY TO COMPUTE, OR CAN BE STORED 

• IRS REQUIRES SCALAR TRIDIAGONAL INVERSION ALONG EACH 
GRID LINE 

ADDS APPROXIMATELY 20% TO CPU TIME BUT BOOSTS TIME STEP 
BY 652% FOR A NET DECREASE OF 4G0% IN CPU TIME 

• EXPLICIT R-K SCHEME WITH IRS COMPETETIVE WITH IMPLICIT 
SCHEMES FOR UNSTEADY PROBLEMS 


CD- 89-39444 


SSME STATOR AND ROTOR GRIDS 


CD-89-39447 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


DETAIL OF STATOR/ROTOR GRID OVERLAP 



UNSTEADY LIFT COEFFICIENTS ON 
STATOR AND ROTOR 


STATOR 



ROTOR 


NUMBER OF PITCH ROTATIONS 
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PRESSURE DISTRIBUTION ENVELOPE ON STATOR 


PS/PO 

STATOR 



PRESSURE DISTRIBUTION ENVELOPE ON ROTOR 


PS/PO 

ROTOR 
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ABSOLUTE MACH NUMBER CONTOURS 



CD-89-37452 


ENTROPY CONTOURS 



CD-09-39453 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


VELOCITY VECTORS AROUND ROTOR 



CD-8?-3?434 


UNSTEADY ROTOR/STATOR INTERACTION CODE 

CONCLUDING REMARKS 

. TIME- ACCURATE ROTOR/STATOR INTERACTION CODE HAS REEN DEVE 
OPED 

. CODE IS APPLICABLE TO AXIAL OR RADIAL MACHINES 


. VISCOUS RESULTS FOR SSME TURBINE SHOW 
UNSTEADY LOADING k SEPARATION, k MIGRATION OF STATOR WAKES 


. IMPLICIT RESIDUAL SMOOTHING 
DOES NOT EFFECT FORMAL TIME ACCURACY OF R-K SCHEME 

ADDS 20% TO CPU TIME 

INCREASED TIME STEP BY A FACTOR OF 6.0 

RESULTS IN A NET DECREASE OF ABOUT 450% IN CPU TIME 

MAKES EXPLICIT R-K SCHEME COMPETETIVE WITH IMPLICIT SCHEMES 
FOR UNSTEADY PROBLEMS 
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